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Abstract 

Patterns of nestedness and specialization asymmetry, where specialist species interact mainly 
with generalists while generalists interact with both generalists and specialists, are often observed 
in mutualistic and antagonistic bi-partite ecological networks. These have been explained in 
terms of the relative abundance of species, using a null model that assigns links in proportion 
to abundance, but doubts have been raised as to whether this offers a complete explanation. In 
particular, host-parasite networks offer a variety of examples in which the reverse patterns are 
observed. 

We propose that the link between specificity and species-richness may also be driven by 
the co-evolution of hosts and parasites, as hosts allocate resources to optimize defence against 
parasites, and parasites to optimize attack on hosts. In this hypothesis, species interactions are 
a result of resource allocations. This novel concept, linking together many different arguments 
for network structures, is introduced through the adaptive dynamics of a simple ecological toy 
system of two hosts and two parasites. 

We analyse the toy model and its functionality, demonstrating that co-evolution leads to 
specialization asymmetry in networks with closely related parasites or fast host mutation rates, 
but not in networks with more distantly related species. Having constructed the toy model and 
tested its applicability, our model can now be expanded to the full problem of a larger system. 
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1. Introduction 



What drives the association between the host-specificity of parasites and the parasite species- 
richness of the hosts that they infect? Generalist parasites are often found in hosts with both 
high and low parasite species diversity, while specialist parasites are found mainly in hosts with 
a rich diversity of parasites (Poulin, 1997; Vazquez et al. 2005). This association, between 
the specificity of parasites and the species-richness of the hosts that they infect, is known as 
specialization asymmetry (Vazquez et al. 2005), or in a slightly stronger case as nestedness, 
when those species occurring in a species-poor assemblage form a non-random subset of those 
assemblages with a higher species richness (Poulin and Guegan 2000). These patterns can be 
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seen from parasites in fish species (Poulin, 1997) to fleas and their hosts (Vazquez et al. 2005), 
although there is some debate on the extent to which this is evident, and many counterexamples 
exist (IPoulinl [19971 [2007] IValtonen et aU|200lb. 



Nestedness is also evident in many other food webs (Ings et al. 2009), and there is mixed 
evidence for whether the addition of parasites to these webs should increase or decrease their 
relative nestedness (Hernandez and Sukhdeo, 2008; Lafferty et al. 2006). Nestedness is a feature 



which is particularly prevalent in mutualistic networks (Bascompte et al. 2003; Vazquez and 



Aizen , 2003 , 2004 ) , an interesting fact given the complete reversal in interaction types involved 



when compared to host-parasite networks. Nevertheless, in order to fit parasites into food web 
models, it is important to understand the forces behind such structural traits as nestedness, 
particularly as antagonistic networks are instead generally expected to be compartmentalized 
( |Bascompte"l|2010| |Thebault and Fontainel|2008[ |Thompson||2005[ but see |Flores et aE1|2011| for 
the effects of scale-dependence). 

What is clear is that structural patterns of species interactions in ecological networks, such 
as nestedness and anti-nestedness, are not random (Joppa et al. 2010). Vazquez et al. (2005) 
constructed a null model based on host abundance to account for such patterns, yet there may 
be many more explanations for the link between the specificity of a parasite and the parasite 



species-richness of the hosts that it infects (Lewinsohn et al., 2006, Poulin, 2007). This does 



not, for example, explain the frequent occurrence of anti-nestedness. Further explanations for 



nestedness in networks include complementarity (Rezende et al. 2007), based on phenotypic 



matching between species, and competitive load (Bastolla et al. 2009), based on a new species 
entering a network targeting a host with less competition provided by resident parasite species. 

Another possible driving force behind the link between specificity and species richness is 
related to the levels of defence that a host exhibits. For example, avian fleas with different levels 



of specificity target hosts with different levels of T-cell mediated immune response (M0ller et al. 



2005). Generalist parasites target hosts with weak levels of immune response, while parasites 



with fewer host species exploit those with both strong and weak immune responses. This also 
has an effect on the parasite species-richness of host species, with hosts with stronger immune 
responses being parasitized by a greater number of species. In this instance, then, the specificity 
of parasites and specialization asymmetry appears to be related to the host immune response 
(M0ller et al. 2005). Here, this idea is turned around slightly, but the concept of a relationship 
between host response and specificity is maintained. 

It is becoming increasingly apparent that the co-evolution of species is an important driving 
force in host-parasite relationships (Best et al. 2009). It is also widely acknowledged that trade- 
offs in resource allocation are responsible for much evolutionary drive in parasites, such as the 
link between virulence and transmission) May and Anderson, 1983). In our paper it is postulated 
that the link between specificity and richness may be driven by the co-evolution of both hosts 
and parasites when balancing the allocation of resources (see Poulin and Morand 2004). These 
resources are devoted in different degrees to interactions with one species versus another, either 
for infection or for preventing infection. 

As an example, the influenza virus binds to cell-surface oligosaccharides via a sialic acid recep- 
tor. The receptor type may have one of two conformations: Neu 5Aca(2,3)-Gal or Neu5Aca(2,6)- 
Gal. A host species may have one of the two or both (Cobey et al. 2010), but the virus must 
adapt to one linkage type at the expense of the other. Other examples include phenotypic trait 
matching, such as the shape of mouthparts of ectoparasites ( Graham et al. , 2009 ) or nectar 
holder size and shape in mutualistic pollinator networks (Vazquez et al. 2009). From the host 
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Figure 1: Possible interaction matrices for a 2-host (y-axes) 2-parasite (x-axes) network. A black 
square represents an interaction between two species and a white no interaction. Plots demonstrate 
specialization asymmetry: (a) a generalist and a specialist parasite and (b) two specialist parasites 
in the same host, and specialization symmetry: (c) two specialist parasites in opposite hosts and (d) 
two generalist parasites. 



perspective, a behavioural example involves reindeer and other herd animals, which often group 
in a reaction to parasitism by biting flies, reducing their chances of attack. This does, how- 
ever, lead to an increased risk of exposure to other pathogens which rely on host density for 



transmission (Hart, 1988). 



Using the above concepts, the following is proposed: if an infection is more prevalent, then 
a host will have a higher likelihood of coming into contact with it and adapting to fight it, 
allocating a greater amount of resource to this and increasing its immune response. A host has, 



however, limited resources with which to do so (Poitrineau et ah, 2003; Poulin and Morand 



2004). Similarly, a parasite may infect multiple hosts. It will, however, be better adapted to 



infect some than others, and again there will be a trade-off in terms of its efficiency in infecting 
a host species (see Poulin|1998 ) . It is therefore assumed that both host and parasite species trade 
off their resources between those species that they target. This trade-off aims to incorporate all 
of the ideas discussed above. A more abundant species will provide more available hosts for a 
parasite. A lower competitive load will encourage infection of that host as an untapped resource. 
Lastly, complementarity will ensure that a species will target another with complimentary trait 
values (parasites will infect hosts which are more vulnerable to them, while hosts will protect 
themselves against parasites which are more of a threat). 

With this in mind, we create a toy model here which investigates the co-evolution of trade-offs 
in a dynamical manner, for a four-species system containing two species of each type; hosts and 
parasites. This may be thought of as a cluster of species forming a 'compartment' in a larger food 



web, and hence more general results may be inferred from the results obtained (see Joppa et al. 
2009). Although this is simply a toy problem, figure [T] gives an indication of what specialization 
asymmetry might look like in this case, and hence the patterns we might expect the model to 
show for nestedness and anti-nestedness in a larger system. It is important to stress that in order 
to fully understand the influence of these trade-offs on nestedness, the investigation of a larger 
system is necessary (in prep.). 

We outline a model describing the system and its mathematical motivation, following this with 
a brief analysis and an investigation of the co-evolutionary dynamics. The results are compared 



to a previous null model based on abundance (Vazquez et ah, 2005), in order to demonstrate 
the compatibility of this theory with our model. 



2. Model 

A standard susceptible-infected system is assumed, with two species of hosts and two of 
parasites. This model has the potential for both specialist and generalist parasites (in one or both 
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hosts respectively) as well as species-poor and species-rich hosts (containing neither parasite, 
one only or both). Si refers to susceptibles of host type i, while Ijj refers to infecteds of host 
type i with parasite type j, where in this instance i,j £ N 2 . We then have 

dS 

=aiNi - aijCijSiFj - UiSiNi, 
jen 2 

=aijCijSiFj - njlij - coilijNi. (2.1) 

The model contains birth (aj) and death (wj) rates dependent on the host species i, as well as 
infection- related death; death rate 7^ of host species i due to parasite species j. Ni = Si+In+Ii 2 
represents the total population size of host species i. 

The maximum force of infection Fj of parasite species j on host species i is given by 

feeN 2 

where /3*- fc is the pairwise potential infectious contact rate for the transfer of parasite j from host 
k to host i. In our model, the actual force of infection is assumed to be moderated by the 
strategies adopted by the parasite j and the host i. It is given by 

G) = OijdjFj = a ijCij ((3fl l3 + /3fl 2j ), (2.2) 

with < djj < 1, < Cij < 1. Here Oj,- is a parasite-related trait defining the relative probability 
of success of parasite j's attack on host i, and Cjj is a host-related trait defining the relative 
probability of failure of host i's defence against parasite j. All else being equal, parasites benefit 
from values of <Zy that are as high as possible, while hosts benefit from values of Cjj that are as 
low as possible. 

We assume, however, that each parasite species j has a fixed amount of resource to allocate to 
infection, and that there is therefore a trade-off between the strength a\j of its attack against host 
1 and the strength a 2 j of its attack against host 2. This trade-off is assumed to be a decreasing 
function, which is species specific and is not dependent upon the population or environment. A 
host species i, meanwhile, varies strategy Cy in order to reduce transmission of parasite j, and 
a similar trade-off is presumed. 

Transmission of infection to a susceptible host then depends on the actual force of infection, a 
measure of both the propensity of the parasite to infect that host, as well as the host's propensity 
to defend itself against the parasite. 



In reality, it is difficult to determine the shapes such trade-offs take (Best et al. 2009). A 



general trade-off shape therefore allows for a greater understanding of different possible evo- 



lutionary outcomes (Kisdi, 2006). Points at which an evolutionary stable strategy (ESS) may 



exhibit evolutionary branching also often depend on the nature of the trade-off function pre- 



sumed; more specifically, whether it is concave or convex, and the extent of this (Kisdi, 2006). 



Levins' fitness set approach (Rueffler et al. 2006) has previously influenced intuitive thoughts 



on the effects of these shapes under normal evolutionary conditions, with a generalist expected 
if a trade-off is convex, or weak, and either specialist expected if a trade-off is concave, or strong. 

The trade-off shape is determined here by a species-specific power (9j for parasite species j 
and (pi for host species i), where 
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a 2j = (1 - a£) e J and c i2 = (1 - cf{) * 
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Figure 2: Examples of trade-off shapes for trait values of parasite species 1, where the trade-off is 
either strong (6*i = 0.5) or weak (9i = 2) 



For parasite species j, 9j < 1 implies a strong trade-off, and 6j > 1 implies a weak trade-off 

(see figure [2]). A parasite is a perfect generalist if a\j = a 2 j = (0.5) e i , which is henceforth termed 
the neutral point. It is a complete specialist if ay = and a 2 j = 1, or vice versa. This is similar 



for host species. Note that the tendency is for trade-offs to be strong, not weak (Ruemer et al. 



2006). 



The second derivative of the trade-off function, the degree to which it is concave or convex, 



can have an enormous effect on the possibility of evolutionary branching (see Kisdi 2006 ) . The 
sign of this depends wholly upon the value of the species-specific trade-off shape, and so the 
system will exhibit very different behaviour around this point. The more concave a function is, 
the more likely the system is to exhibit branching (Kisdi 2006). Note that the second derivative 
is maximized (either positive or negative depending on the trade-off shape) for parasite j at 

a±j = (0.5) ^ , where 



d 2 a 2 j 

~daj: 



(0j ~ 1), 



and similarly for host species. This affects the direction in which trait-values mutate, which will 
be discussed in greater detail in the following sections. 



3. Analysis 



This model is analysed using adaptive dynamics. This assumes rare mutants with marginally 
different phenotypic effects to residents, which may then invade the population if their growth 



rate is positive in an equilibrium environment (Dercole et al. 2003). This approach does as- 



sume clonal reproduction, but the results still hold for random mating in monomorphic diploid 



populations with polygenic traits if mutations are rare with small phenotypic effect (Law et al 



2001 Ruefner et al. 2006). The growth rate of a mutant phenotype while rare is termed the 



invasion fitness, and is important in determining whether that mutant may invade, and poten- 



tially replace, a resident (de Mazancourt and Dieckmann 2004). The population evolves in the 



direction of the fitness gradient (the change in the invasion fitness with respect to change in 
the mutant trait value) as successive mutations occur and then spread through the population 



(de Mazancourt and Dieckmann, 2004). The invasion fitness can then be used to discover singu- 



larities, where the fitness gradient of any local mutant is zero, and to investigate the nature of 



these singularities (Geritz et al. 1998). An example of the derivation of the invasion condition 
is given below for parasite 1 in a system without co-infection or recovery; the method is similar 
for the second parasite species. 

3.1. Parasite invasion conditions 

In determining the invasion fitness of a mutant parasite in the population, here a mutant 
I[i with strategy (a'ii,a 21 ) of parasite 1, with the resident population at a stable, non-trivial 
equilibrium (Si = S*, Iij = I?), the linearized dynamics of the mutant are expressed as follows: 



dt 
dl^ 

dt 



a' 11 cnSlF 1 - 7ii-fn - uil'nNii 
a / 2i c 2i5 , 2 i? i' - 721^21 - ^2/21^2*- ( 3 - 3 ) 



In the full system of equations for the dynamics, including the resident dynamics, the resident 
is at a stable steady equilibrium. The Jacobian of the system can be therefore split into the 
original dynamics and a separate submatrix determined by the mutant, given by A\. 

Defining 8\\ = a' u — an and S21 = a' 21 — 021, both small, we get 



(3.4) 



Ai(l,l) Ai(l,2) 
Ai(2,l) Ai(2,2) 

where the equilibrium of the system is used to guarantee that 

I* 1 ( lil + L H N*) = a il c il S?F?, 

and hence 

Ai (n,n)I n i = 5 n ic n iS n f3i' 1 I n i — a n ic n iS n (3™ n 1 / ra _ 11 , 

Ai(n,n- 1) = (<5„i + a nl )c nl S* n ^ n ~ l ■ (3-5) 

The stability of Ai is then used to indicate the potential for invasion. As it is difficult to 
interpret anything useful from the eigenvalues, a sign-equivalent proxy can be found for the 
growth rate of the mutant by investigating the trace and determinant of A\. The determinant 
is given by 

IiiiJi det(Ai) = -ciiC2iS , *5 , 2-ui((5ii,aii,52i,a2i), 
where, to leading order of 5, 

ui(6u,axi,52i,a2i) =f3\ 2 f3j 1 /n-^i^ii^i +011^21) 
+ Pl 1 p\ 1 {r il ) 2 5 ll a 21 

+ P{ 2 Pl 2 {r 2l ) 2 a ll 5 21 . (3.6) 

When a'u = an, then tr(Ai) < and det(Ai) = 0, as 8 n \ = 0. For the small perturbations 
resulting from mutation we can therefore rely on the determinant condition for stability analysis, 
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as the trace will remain negative. The cases for I*i = and I|i = are discussed later. The 
determinant will be negative, and hence the system is unstable and the mutant invades, if 

ui(£n, a n , 621, 021) > (3.7) 

Interpreting this condition, invasion of a mutant parasite can be seen to be driven by the 
following; the first term is due to inter-species transmission, and is minimized at the neutral point 
an = OL21, hence driving towards generalism. This is due to an increased invasion probability 
when a mutant is investing more equally than the resident, as (Sna 2 i + 011^21) is larger. When 
011 = fl 2i then this term is zero, and no possible mutant can invade due to this. When an < 

(0.5) e i , this implies that 021 > (0.5) e i . Hence we require a mutant with a larger value of an to 

j_ 

promote invasion, and vice versa for an > (0.5) "i . So this term stabilizes at and attracts to the 
neutral point, promoting generalism. 

The final two terms compare the use of the two available hosts. For example, if transmission 
due to host 1 (/3 2 P\ (-^Ti) 2 ) i s high, then a mutant with larger a n will invade, ensuring host 
1 is utilized. These terms demonstrate an increased invasion potential if the mutant increases 
infection of the species with higher infection rates (/S^ 1 /^ 1 compared to Pl 2 f3 22 ) and on which 
the resident relies more (whether I n is greater or less than I 2 i), as <5n and 621 have opposite 
signs. This ensures that the mutant is spread as much as possible to susceptible hosts. 

Note that the case where either I u = or 7|i = sees the parasite shy away from a completely 
protected host. These points are only obtainable, for a four species system where the maximum 
force of infection is not zero, if the actual force of infection is zero, i.e. aycii = or 021021 = 
as appropriate. For the case where I u = 0, the submatrix of the Jacobian is given by 

-7ii-wu(Si +^12) 
(<$2i + a 21 )c 21 S* 2 (3 21 5 2l c 2l S* 2 Pf 

The mutant can thus invade only if 5 2 \ > 0, i.e. if a' n < an. A similar situation arises for 
I 21 = 0. If both I n and I 2 \ = then the system is at a trivial equilibrium, in which case no 
mutant can invade. 

3.2. Host invasion conditions 

Similarly to the parasite case above, a mutant population (S[, with trait value c^- of host 
1 is introduced at low densities to the resident population at equilibrium (Si = S*, Iij = 
The dynamics of the subsystem are given by 

f]QI 

= ai N[ - a u c' n S[Fl* - a l2 c' l2 S' l F 2 u - ^[N?, 

dV 

=an C ' 11 S 1 F 1 1 * - 711/n - u x l' n Nl 

dV 

=ai2C , i2S 1 F 2 1 * - 712/12 " wili 2 JVj. (3.9) 



The equilibrium conditions are taken from equation 2.1, and we define en = c u — en and 
ei2 = c'12 — C12, both small. Again we investigate the submatrix of the Jacobian, given by 

/ Cx(l,l) ai qi \ 

Ci= a n (en + cn)Ft - a ^ F '* . 

v ai 2 (e 1 2 + ci2)F 2 1 * _™pi! 



where 



/ at* \ 

Cx(l, 1) = olx [\ - -^j - a n e u Ft - a 12 e 12 F 2 u . 

Now, as the subsystem is three-dimensional, the eigenvalues are given as solutions to an 
equation of the form A 3 + 61 A 2 + 62 A + 63, where 

61 =-tr(Ci), 
2b 2 =tr 2 (d)-tr(C?), 
63 = - det(d). 

For stability, it is required that 61,62,63 > and 6162 > 63. When c' u = en, both 61 and 62 
are positive and 0(1), as e = 0. On the other hand, 

63 = ~ T J~ T * anauSlF^* F 2 u wi(en, c n , e u , cu) , 
-'n-'i2 

where 

wi(eii,cn, ei2, C12) =ci2enai/ n + c\\t\ 2 a\l\ 2 

- cuc 12 auSlFl* en 

- c n ci 2 a 12 SlF2*e 12 . (3.11) 

Now this is 0(e), and is zero when c' n = en, hence 6162 > 63. The criteria for instability (and 
invasion) when mutations are small can therefore be reduced to 

u>l(eu, en, ei2, ci 2 ) > 0. (3.12) 
Interpreting this condition, it can be seen that mutant invasion depends on a balance of 



terms. The first two terms from equation 3.11 push the host towards a generalized defence. For 
example, a particularly large value of c\\ makes the second term likely to be larger, and hence 
a mutant with a smaller trait value c' n will invade, to ensure that e\ 2 > 0. 

The final two terms decrease in importance {cuc\ 2 decreases), compared to the initial terms, 
as the host specializes. If the pressure from parasite 1 is higher (the actual force of infection is 
higher) then the first of these terms will be larger. A host mutant with a smaller trait value will 
therefore invade, in order to make en negative and the sum of the final two terms positive. In 
this way the mutant host protects itself against parasite 1 to a greater extent. These terms then 
account for the pressure that a parasite places on the host, and the host's reaction to this. 

The case where either I n or I\ 2 is zero, which again occurs if a\\c\\ = or a\ 2 c\ 2 = as 
appropriate, shows that a host will not defend itself against a non-threatening parasite, and if 
the system is at a trivial equilibrium, then no mutant may invade. For example, if / n = 0, then 



equation |3.10| is given by 
Ci 



I Ci(l,l) ai ai 

-711 - 

v a X2 d l2 Fl* 



-711-wiJV? | ; (3.13) 

1* n ai2Ci 2 S*Fi* 



where C\(l, 1) = aj ( 1 



a i2 e i2-p2*- I n this case, a stronger than necessary condition is 



that a mutant will always invade if 012612-^2* < 0> i- e - ^ c ii > c n- A similar situation arises for 

T* — 
i 12 — u. 

The conditions for equilibrium can be independently verified using the next-generation tools 



outlined in Hurford et al. (2010). For p(M) defined to be the spectral radius of matrix M, and 
s(M) the spectral bound, the Jacobian C\ is decomposed such that 



(3.14) 



for F\ and V\ satisfying s(— Vi) < 0, V^ 1 > and F\ > 0. i^i and V\ are taken as matrices 
representing the appearance and disappearance of individuals from the system respectively. The 
invasion condition is then given by 



p(F 1 V^ 1 ) > 1, 



For further details see Hurford et al. (2010). In this case, we take 



and 



Vi 



\ 



a\ a,\ a\ 





Vi(l,l) 
-a u (eii + c n )Fl* 

-ai2(ei2 + cn)F%* 



(3.15) 



(3.16) 













(3.17) 



/ 



Nl 



012^12-^2*- After some straightforward algebra, this yields 



where Vi(l, 1) 

an identical condition to equations 3.11 and 3.12 This approach is not useful for the parasite 



invasion conditions, as individuals enter the system through more than one class (see Hurford 
leTaT1l20T0l 

From the initial analysis, therefore, it appears that parasites infect more vulnerable hosts, 
but also aim to be generalists. Hosts trade-off between the pressures exerted by different parasite 
species, aiming to lower this as much as possible, but similarly aim to generalize their defence. 
Note that the cases for mutants of parasite and host species 2 are similar to the above. This is 
all as expected, and demonstrates the biological validity of our trade-off. 

Simplifications of the above scenario demonstrate that the co-evolution of all four species 
is vitally important, as expected. If one or more species is missing from the system then the 
trade-off ensures that the dynamics are trivial, with those that have two antagonists splitting 
their resources as before, while those that have only one concentrate their resources on that 
one. This demonstrates the importance of such a system, incorporating both multiple hosts and 
parasites, unlike many previous approaches taken when modeling host-parasite systems. 



4. Results 



Simulation of the results follows the method of Dieckmann and Law (1996), concerning the 
frequency and impact of selection, where evolutionary dynamics occur at a much slower rate than 
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(a) 



(b) 



Figure 3: Pairwise invasion plots for parasite 1 in a symmetric environment, with intra-species 
transmission higher than inter-species transmission. Mutant trait values are on the y-axes, while 
resident trait values are on the x-axes. Trade-off shapes are (a) strong for the parasites and (b) weak 
for the parasites, while linear for the hosts. 



population dynamics (see Drossel and McKane|2005 This relies on the derivative with respect to 
the mutant trait- value of the growth rate of the mutant in a population of residents, given by the 



dominant eigenvalue for each species from matrices 3.4 and 3.10 As a result of the appropriate 



eigenvalue equations, discussed in sections |3.1| and |3.2[ and noting that A must be small, the 
dominant eigenvalues may be approximated up to some positive multiplicative coefficients by 
— det(Aj) and det(Cj) respectively. We therefore use the following general equation for species 
k with trait- value Sk and the appropriate approximation for the eigenvalue E^: 



d_ 
df 



Q 
OSi„ 



(4.18) 



where fk is the evolutionary rate coefficient (see Dieckmann and Law 1996 This is calculated 



for an environment determined by resident trait- values s for all species. 

Changes to trait values with time are then investigated, where species can co-evolve. It is 
presumed that (pi = 1 for hosts, so the trade-off shape is linear, and initial trait values are at the 
neutral point (0.5 here), in order to attempt to separate the effects of hosts and parasites. There 
are a number of possible cases for different 9 values; trade-off curves of both parasite species 
may be either concave or convex [9 < 1 or 9 > 1 respectively). 

A look at the pairwise invasion plots for cases with different trade-off strengths demonstrates 
the outcomes that we expect (see figure [3]). Pairwise invasion plots (PIPs) indicate when a 
mutant can invade (shaded) or not (white) depending on its trait value relative to the resident 
trait value. This property depends on the mutant's invasion fitness at low frequency in a resident 
population. The population evolves as small mutations occur which move the population off the 



diagonal (Geritz et al. 1998). If the mutation is successful (i.e. the mutant is in a shaded area) 



then the mutant population grows and it displaces the resident, to become the new resident 



(Geritz et al. 1998). In figure |3(a) here, for example, this occurs until the population's trait 



value reaches an extreme, depending on which side of the central point the resident trait value 



begins. Species with a strong trade-off are likely to evolve to be specialists (figure 3(a)), while 
those with a weak trade-off are expected to evolve towards generalism (figure |3(b) ) . 



Figure [3] demonstrates the evolution of one parasite only, in a static environment where no 
other species evolves. In a full analysis this will not be the case. This motivates us to follow 
the co-evolution of trait values for all four species simultaneously. Results are demonstrated in 
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Evolutionary Time 



Evolutionary Time 



(a) Weak trade-off for parasites (Oj — 2) 



(b) Mixed trade-offs for parasites (9i = 2, 62 = 0.5) 



Figure 4: Examples of the evolution of trait values with time in a symmetric environment, with 
intra-species transmission higher than inter-species transmission. Rates of change of trait values are 
calculated directly from the growth rates of mutant traits in a resident population. Trade-off shapes 
are linear for hosts and (a) weak or (b) mixed for parasites. The inset in (b) shows the initial dynamics 
of the system before the slower host mutations have had an effect. Note that an and 021 represent 
trait values for parasite species 1 in hosts 1 and 2 respectively, while 012 and 022 represent the same 
for parasite species 2. 



figures [4] and [5] 

In figures [4] and [5j although mutation rates are taken to be identical for hosts and parasites, 
the growth rate of mutant hosts in a resident environment is significantly slower than that of 
parasites. The host trait values are not presented here, although their importance is discussed 
later. 

Again we assume a linear trade-off for hosts, which have initial trait values at their respective 
neutral points. In figure [4(a) where #i,#2 > 1 ; generalist parasites evolve. Here the trait values 
of parasites can be seen to evolve towards their neutral points. Note that the system cycles 
around the neutral point until both parasites are perfect generalists. This corresponds to figure 



1(d) 



For the case where the trade-offs for the parasites are strong for one parasite and weak for the 
other (figure |4(b)[ ), we obtain the coexistence of a relative generalist and an extreme specialist 
in the same hosts. Note that this system takes much longer over evolutionary time to equilibrate 
than other cases, due to the slow rate growth rate of mutant host populations. Hence we see 
different dynamics over short (figure |4(b)| inset) and long (figure |4(b)[ ) time-scales. Before the 
hosts are able to react to the presence of the parasites, we see a generalist host that is more 
focused on the opposite host to the specialist (trait values an and 022 are higher), whereas 
once the host mutant populations have had an effect, we see both the relative generalist and 
the specialist parasite are more focused on the same host (trait values an and 012 are higher). 
As a host which contains only a generalist would be able to focus its defensive efforts on that 



parasite, and hence, as can be seen from equation |3.10[ the parasite would not target it, we 
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(b) Parasites in opposite hosts 

Figure 5: Examples of the evolution of trait values with time in a symmetric environment, where 
the trade-off shapes for both parasites are strong (6j — 0.5). Different cndpoints occur as a result of 
the initial equilibrium values for susceptible and infected hosts. These result in parasites in (a) the 
same or (b) opposite hosts. 
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see only comparative generalists here. This system demonstrates specialization asymmetry, and 
indeed nestedness (as far as that is plausible in such a small system) over longer evolutionary 
time-scales, corresponding to figure l(a)| 

For the case where both 9j < 1 (figure [5]), specialists always evolve. These can evolve to be 
in the same hosts (i.e. an and a\2 evolve to the same extreme, figure 5(a)), or in the opposite 
host (an and a\2 at the opposite extremes, figure [5(b)] ) . These scenarios bear resemblances to 
specialization asymmetry (figure 1(b) ) and compartmentalization (figure 1(c) where the network 



is split into separate sub-networks that are not linked to one another) respectively. 



4-1. Initial trait values 

From the parasite conditions for invasion, it can be seen that the behaviour of the cases 
depends heavily on the initial equilibrium conditions, which are a result of the initial trait 
values. The behaviour of each scenario pivots around which side of case-specific points the 
initial trait values lie, similarly to which side of an evolutionary stable strategy an initial point 
lies in a pairwise invasion plot (see figure [3]). Each case has only a limited number of evolutionary 
end points for trait values, and evolves to these. An analysis of the initial trait values shows 
interesting results. 

For the case of a weak trade-off for both parasites (figure |4(a)| ) , although generalists occur 
almost exclusively in "opposite" hosts, the trait values of the parasites evolve to be so similar 
that only a very careful inspection can detect the difference. In other words, two generalists 
occur, each slightly more dependent on a different host. If one parasite has a strong trade-off 
and the second a weak trade-off, a relative specialist and a relative generalist occur respectively, 
in the same host (figure |4(b)[ ),This demonstrate specialization asymmetry, and occurs for all 
initial trait values. For the case where both parasites demonstrate a strong trade-off (figure [5]), 
figure 6(a) shows the initial trait values which lead to specialization asymmetry. 



Increasing the hosts mutation rates (/& from equation 4.18) serves to alter the shape of the 
curves towards that found in figure |6(b) The much faster generation time of parasites has been 
used in the past to justify the study of their evolution alone, as opposed to a full co-evolutionary 
system (Vazquez and Aizen, 2004[ ). Even including co-evolution, however, it has been shown by 
Best et al. (2009) that different mutation rates do have a significant effect. In our simulations, 



where the growth rates of mutant host populations in resident environments appear significantly 
slower than those of mutant parasite populations, increasing the mutation rates of hosts in 
comparison to parasites greatly increases the chances of specialization asymmetry occurring for 
a wider range of initial trait-values when a strong trade-off exists for parasites. This indicates 
not only the circumstances under which we might expect to see asymmetry, but also a possible 
line of experimentation to take in order to validate our results. 

The addition of both recovery terms and co-infected classes to the system, computed numeri- 
cally, demonstrated increased likelihoods of specialization asymmetry occurring, again dependent 
on both the initial trait- values and relative mutation rates of species (unpublished results). 



4-2. Abundance 

In our toy model we included both abundance and phenotypic matching as motivators for 
nestedness. In order to justify the claim that the model is compatible with the concept of 
abundance as a driving force (see Graham et al.|[2009 Poulin and Morand||2004 Vazquez et al. 
2009), we investigate here the manner in which hosts of different abundance influence 
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(a) Low host mutation rate (b) High host mutation rate 

Figure 6: A sketch of the end-points of evolution, in terms of species-richness, for different initial 
trait values for a strong trade-off (Oj = 0.5) for parasite species. Black areas denote those initial 
trait values that evolve to have parasites in opposite hosts, while white areas denote those for which 
parasites evolve to be in the same host for (a) a low host mutation rate and (b) a high mutation rate 
for hosts. Arrows indicate the direction of evolution of trait values. Note that the trait value for each 
parasite with respect to the second host is determined by the trait value for the first, and hence is 
not included in these plots. 
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the model outcomes. This is especially important over evolutionary, as opposed to ecological, 
timescales. 

There are two aspects of abundance to be compared to the model presented here. Firstly, 
does increased abundance of a species indeed lead to a higher number of links occurring, and 
secondly, does asymmetry in the assignment of links lead to nestedness? The first of these is 
compared to the model, while the second is a question that remains independent of the model. 

Firstly we note that, if nestedness is a result of abundance then we would expect it to occur 
in mutualistic and predator-prey webs, as the number of links of a species is associated solely 



with its abundance and not with the nature of those links. This is indeed evident (see Montoya 



et al. 2006), particularly in mutualistic networks, which display more nestedness than would 



be expected from a random, bipartite network (Bascompte et al. 2003). However, if this is the 



case then why would certain web-types display more nestedness than others (see Hernandez and 



Sukhdeo 2008, Lafferty et al. 2006)? This promotes the thinking that abundance alone cannot 



explain patterns of nestedness in ecological networks. A further question concerns whether or 
not abundant species are more likely to have links simply as they are more likely to interact 
with other species due to their abundance, or whether this is due to the inherent benefits of 
interacting with a more abundant species. Our model attempts to address this, proposing that 
there is an evolutionary advantage to interacting with a more abundant species. 

For our trade-off model, note in equations 3.6 that the invasion potential of a parasite de- 
pends on the relative abundances of the different host species. The terms In and I21 are the 
combination of infection prevalence and total population size, and hence reflect that a larger 
population will increase invasion potential. This was also computed numerically, where it was 
observed that, in a symmetric environment, increasing a species' population size could drive a 
parasite to preferentially interact with that species (see figure [7|) . Here we look at parasite 1 
interacting with host 1 and see that, although this is still dependent on initial values for traits, 
increasing the abundance of the host species clearly increases the likelihood of the parasite 
evolving to interact with that species. 



5. Discussion 



This model provides evidence for both specialization symmetry and asymmetry, but it is diffi- 
cult to comment on nestedness from this. A similar model might be sufficient for larger networks, 
where a species-poor host could still contain more than one parasite, allowing parasites to split a 
host's defence. Additionally, if hosts were to trade defensive properties in an alternative manner 
(for example through reduced birth-rate, or the inability to reduce transmission completely to 
zero) this could allow generalists to exist alone in a species poor host. 

There is a great deal of mixed evidence, from a number of different ecosystems under varying 



conditions, for nestedness in host-parasite interactions (Graham et al. 2009). In those systems 



in which nestedness occurs, generalists will be in comparatively species-poor hosts, as they may 
occur in both species-poor and species rich hosts. In this sense the model can be related to 
nestedness, although its size makes a full comparison unreliable. The model also supports the 
idea that a strong trade off leads to specialists, while a weak trade-off promotes generalism, in 



accordance with the Levins' fitness set approach (Levins 1962, cited in Rueffler et al. 2006 



From the results obtained here it is evident that initial trait values are very important to the 
final equilibrium of a population. This may be useful in predicting responses when a species is 
added to or removed from a system, as such an event will be followed by evolution of the system 
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Figure 7: Regions of initial trait values for both parasites which, under co-evolution with a strong 
trade-off, lead to parasite 1 occuring in host 1. Filled areas denote regions of initial trait values which 
see parasite 1 evolving to be in host 1, and succesively lighter shades denoting the increase in size of 
the region as the relative abundance of host species 1 increases. 



in a direction dictated, to a certain extent, by the equilibrium values of the system prior to the 
alteration. 

One crucial factor concerning this model is that it is dynamical. A structural property can 
be described and predicted by this dynamical model, explaining an aspect of static models that 
cannot be explained through a simple trophic hierarchy. This helps in clarifying how parasitic 
associations may be motivated, and, to a certain extent, investigates the effects which parasite 
in a host species have on one another. This toy model can now easily be expanded to include a 
more realistic system containing many more species (in prep.). Analysis of a larger system will 
then enable questions on the nestedness of the system as a result of trade-offs to be answered. 



5. 1 . Future 

The primary purpose of this paper was to lay the foundation for a model which could explain 
patterns of nestedness in ecological networks. In order to do so, this model needs to be repeated 
at a larger scale for a host-parasite network. Given the evidence for nestedness in other systems 



(Bascompte et al. 2003; Vazquez and Aizen, 2003, 2004), the adaption of this model to these 



systems, particularly mutualistic networks, could also be used to corroborate any conclusions 
reached (in prep.). 

Many other factors are considered as possible motivators for the species-richness of parasites 



in hosts (see Feliu et al. 


1997 


Marcogliese 


2002 




Montoya et al. 


2006, 


Morand and Poulin 


1998 


Nunn et al. 


2003 


Poulin and Leung 


2011 


Poulin and Morand 


2004 


), which focus on the 



characteristics of the hosts. It is becoming increasingly clear that the interests of the parasites are 
also important factors in this, and, in fact, both motivating factors are likely to be of importance. 
There have been very few models which investigate the co-evolution of a host-parasite system 



using adaptive dynamics (but see Best et al. 2009; Caval and Ferriere 2010), and these focus 



primarily on the discovery of a co-evolutionary stable strategy. In that sense, every additional 
approach to or analysis of a co-evolutionary system adds to the field of co-evolutionary ecology. 
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With the aim of co-evolution in mind, it has been observed that, according to game the- 
ory, predators may be responsible for the presence of additional prey species through induced 
branching (McGill and Brown, 2007). Is it possible that parasites have such an effect? This could 
presumably only happen if parasites exerted similar levels of pressure to predators, which is un- 
likely (Poulin 2010). Branching in our model could, however, potentially explain nestedness, as 
parasites would be found in similar hosts. This would lead to results similar to figure 5(a), an 
aspect which could be investigated further. 

A further step from here is to investigate the effect that the position of a host in the network 
as a whole has on its parasites, and how this fits in with the observations made here. The 
position of a host species in the network is considered a potential driving factor in determining 
its parasitic composition (Chen et al. 2008 Vazquez et al. 2005). This has been looked at to a 
greater extent than parasite interactions with each other in the past (Cobey et al. 2010), but 
open questions still remain, specifically with larger networks that include both multiple hosts 
and multiple parasites together. 



6. Conclusion 

The results of the model indicate that the hypothesis of resource trade-off driving a link 
between specificity and species richness appears to be plausible. It can certainly be used to 
model interactions between hosts and parasites, which should yield interesting results when 
used on a larger scale. This also highlights the importance of factors such as host mutation rates 
in co-evolutionary systems, even when these rates are low. 

Using such information as our results for the mutation rates and initial trait values, our model 
helps to predict the circumstances under which we might expect patterns such as specialization 
asymmetry to occur. We would predict the presence of specialist parasites in species-rich hosts 
to be more likely if the hosts had higher mutation rates, and in systems in which parasites are 
more closely related, are more likely to originate in similar hosts or appear as generalists. Given 
the relationship between specialization asymmetry and nestedness, we would therefore expect 
nestedness under similar circumstances, and anti-nestedness the remainder of the time. 

This model demonstrates that dynamic co-evolution of the network is vitally important in 
accounting for parasites, as it demonstrates how the dynamics could influence structural proper- 
ties. In particular, it demonstrates the importance of the co-evolution of both hosts and parasites 
in such a scenario. Parasites, therefore, are not a characteristic to simply be transposed onto 
a system with no regard to their effects on one- another. Much like interactions in conventional 
food webs, the influences of different parasites can alter the entire structure of a host-parasite 
network. 
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